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Abstract 

^ ■ We study stabihty of superflow of Bose gases in optical lattices by analyzing the Bose-Hubbard 



model within the Gutzwiller mean-field approximation. We calculate the excitation spectra of the 
homogeneous Bose-Hubbard model at unit filling to determine the critical momenta for the Landau 



1-^ ' and dynamical instabilities. These two critical momenta are shown to approach each other when 

a ' 

Q . the on-site interaction increases towards the Mott transition point. In order to make a direct con- 

O 

nection with realistic experiments, we next take into account a parabolic trapping potential and 

t^ compute the real-time dynamics of dipole oscillations induced by suddenly displacing the trap cen- 

Q«^ ■ ter. We consider the following two cases: standard softcore bosons, whose interparticle interactions 

Tij- , include the on-site one only, and hardcore bosons with long-range dipole-dipole interactions. For 

^— V \ both cases, we show that the dipole oscillation is significantly damped when the maximum local 

momentum exceeds a certain threshold, which quantitatively agrees with the critical momentum 
for the dynamical instability in the homogeneous system. In the case of dipolar hardcore bosons, 
the dynamical instability of dipole oscillations leads to the formation of checkerboard density waves 
in the superfluid phase near the boundary to the supersolid phase. 
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I. INTRODUCTION 



Since the observation of the quantum phase transition from superfluid (SF) to Mott insu- 
lator (MI) [l|, systems of ultracold gases confined in optical lattices have provided an ideal 
venue for studying strong correlation physics |2|. Many interesting properties of strongly 
correlated ultracold gases have been revealed in recent experiments, including excitation 






spectra |3|-l6|, quantum criticality [7|, and particle transport |8l-ll0|. 

In particular, transport properties have attracted much interest, in connection with the SF 



critical momentum above which superflow breaks down. In experiments with ultraco 
in optical lattices, transport is investigated by using a moving optical lattice 9|, 
suddenly displacing the parabolic trapping potential to induce a dipole oscillation 
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In the weakly interacting regime, where the Gross- P it aevskii (GP) mean-field approximation 
is valid, it has been shown that the critical momentum in a trapped system at low tempera- 
tures agrees quantitatively with that for the dynamical instability (DI) in the homogeneous 
lattice system llMlSj , regardless of whether one uses the trap displacement or the moving 
optical lattice. In the intermediate and strongly-interacting regimes, while the agreement 
in the critical momenta of the tr app ed and untrapped systems has been obtained when the 



moving optical lattice is used 
the trap displacement. 
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17(1, it remains unclear whether this is also the case for 



In addition to optical lattices, new possibilities for the study of strong correlation physics 
have been opened up b y th e creation of ultracold atomic gases with strong magnetic dipole- 
dipole interactions [18|, |l9| and gases of polar molecules 20|, |2l|. It has been predicted that 
when a dipolar Bose gas is loaded on an optical lattice, there exist supersolid (SS) phases that 
possess both diagonal (crystalline) and off-diagonal (superfluid) long-range order |22|-|29|. 
Danshita and Yamamoto have studied the critical momenta of dipolar Bose gases in a two- 
dimensional (2D) optical lattices and suggested that some properties of the critical momenta 
can be used to identify supersolid phases 28l|. More specifically, the critical momenta of 
the SS phases are finite in contrast to the insulating phases, such as MI and density wave 
phases, and significantly smaller than that for the SF phase. Since the analyses in Ref. 28 1 
have been done in a homogeneous lattice system, it is important to investigate the critical 
momenta of dipolar Bose gases in the presence of a parabolic trap in order to make the 
suggestion more convincing. 



In this paper, using the Gutzwiller approximation, we study superfluid transport of Bose 
gases in optical lattices with and without the dipole-dipole interactions. First, we analyze 
the critical momenta of Bose gases without the dipole-dipole interactions at unit filling in 
homogeneous ID, 2D, and 3D optical lattices by calculating the excitation spectra. We locate 
the boundaries to the Landau instability (LI) and the DI. Secondly, solving the equation of 
motion numerically, we simulate the dynamics of dipole oscillations of Bose gases confined in 
a parabolic potential both in the absence and the presence of the dipole-dipole interactions. 
For both cases, we find significant damping of the dipole oscillations when the maximum 
local momentum exceeds a certain critical value, and that the critical value coincides with the 
critical momentum for the DI in the homogeneous lattice systems. We also find a parameter 
region in which the dipole-oscillation mode is resonantly coupled with the breathing mode. 
When the dipole-dipole interactions are present and the system is in the SF state close to 
the SS phase, we show that the DI of a dipole oscillation is followed by the formation of 



checkerboard density waves, as predicted in Ref. [5 

This paper is organized as follows. In Sec. [TTl we introduce the Bose- Hubbard Hamiltonian 
and our formulation of the problem based on the Gutzwiller mean-field approximation. In 
Sec. mil we calculate the critical momenta at unit filling in homogeneous ID, 2D, and 3D 
optical lattices from the excitation spectra. In Sec. IIV| we calculate the dipole oscillations to 
determine the critical momenta in the presence of a parabolic trapping potential. In Sec. |Vl 
we analyze the dipole oscillations of dipolar hardcore bosons. In Sec. IVIl we summarize our 
results. 

II. MODEL AND FORMULATION 

A. Bose-Hubbard Model 

We consider a system of ultracold Bose gases confined in deep optical lattices combined 



with a parabolic trap. This system is well described by the Bose-Hubbard model 30 



3l| 



d ^ 



H = -J5^5^(ajaj+e. + H.c.) + -^nj(rij - 1) + ^^(ej - /i)nj, (1) 



j a=l 



where ej = r2|aj — ro(t)p. a and d represent the lattice constant and the spatial dimension. 
The vectors j and e^ denote the site index and a unit vector in the direction a, where 



the directions a = 1,2, and 3 mean the x, y, and z directions, cij (a-) is the annihilation 
(creation) operator and n^ = a-d^ is the number operator at the site j. J, U, fi, and /i are the 
hopping amphtude, the on-site interaction parameter, the curvature of the parabohc trap, 
and the chemical potential, respectively. We fix the trap curvature to be Q/J = 0.01 when 
we consider dipole oscillations of Bose gases in optical lattices in the presence of a parabolic 
trap. We assume that dipole oscillations are induced by displacing the trap center in the x 
direction as t'o(^) = Da6{t)ei , where 6{t) is the step function. We control the momentum 
of Bose gases by changing the displacement Da. Henceforth, we adopt the units in which 
J = h = a = 1 except in the figures and their captions. 

B. Gutzwiller approximation 



In order to calculate the ground states, the excitation spectra, and the real-time dynamics 
of Bose gases in optical lattices, we use the Gutzwiller approximation, in which the many- 
body wave- function is assumed to be a single product of local states as |32l |. 



where |?t,)j represents the local Fock state at the site j and the normalization condition for 
coefficient /j,„ is J2n \fj,n\'^ — 1- Minimizing the effective action, dt ( \1/g 



with respect to /=*„, one can derive the equation of motion for /j „, 



dt 
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-n{n — 1) + (cj — yu)n 



(3) 



where $j = {^g\0'j\^g) = Tin V^/jtn-i/j," i^ the superfluid order parameter. While Eq. ([3]) 
obviously describes the real-time dynamics, it also allows us to calculate the ground states 
by solving it in imaginary-time (33|. A steady solution takes the form of /j,n(^) = /j,ne~*'^-'*, 
where f\n is time-independent and the phase factor Cj\ is given by 



UJ\ 
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/^j 



n|/j,nP. (4) 



Here, /ij = /i — ej is the effective chemical potential at the site j and $j = Tn V^fln-ifi,n 
is the superfluid order parameter for the steady state. We will explain how to calculate the 
excitation spectra from Eq. ([3]) in the next subsection. 



The Gutzwiller approximation has been often apphed to the Bose-Hubbard model in 
order to analyze various phenomena and properties of Bose gases in optical lattices. This 
approximation is more accurate in higher dimensions because of its mean-field nature. In- 
deed, thorough comparisons with experiments have shown that it can quantitatively describe 
several interesting properties in 3D, such as the excitation spectra |6| and the critical mo- 
menta |9|. The qualitative validity of this approximation is thought to hold also in 2D, given 
that it correctly captures the order of the SF-to-MI quantum phase transition. In contrast, 
it is widely known that the Gutzwiller approximation miserably fails in ID except for very 
weak interactions f/ ^ 1. Nevertheless, in the present paper we analyze ID systems using 
the Gutzwiller approximation because the simplicity of ID systems is very useful to illustrate 
basic properties of the superfiuid transport that holds in the level of the approximation. 

C. Excitation spectrum 

Previous works have described prescriptions to calculate the excitation spectra of the ho- 
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34| . We analyze 



mogeneous Bose-Hubbard model within the Gutzwiller approximation 
the excitation spectra not only to investigate the stability of Bose gases in homogeneous sys- 
tems, but also to obtain collective-mode frequencies in trapped systems. While Kovrizhin et 
al. have used an extended version of the prescriptions to calculate the excitations in trapped 
systems 



35| . they have not explicitly explained their formulations. In this subsection, we 
present a detailed explanation of the extended version that allows us to deal with spatially 
inhomogeneous systems. 

We consider a small fluctuation Sf-^^n around a steady solution /j^^i, 

hn{t)= [/j,n + '5/j,n(t)]e-^^J*, (5) 

We substitute Eq. (|5]) into Eq. ([3]) and linearize the equation of motion with respect to the 
fluctuation (5/j ,1. The Bogoliubov transformation on the fluctuation, 

5h,n = «j,ne-^-* + vl^e^-*\ (6) 



leads the following linear equations, 
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-X15Z [{v^^^/j,n-l/j*_e„,„^_l + ^{n + l)(m + l)/j,n+l^*_e,,m+l J "j-e„,m 
m a=l 

+ |y^^/j,n-l/j+e<„m-l + V(n + l)(m + l)/j,n+l^+e„_m+l } «j+e.,m 
+ I Vn(m + l)/j,n-l/j-e,,m+l + \J{n + l)m/j,„+i/j_e„m-l} ^^j-e„,m 
+ { A/n(m^ + l)/j,„_i/j+e„,m+l + \/{n + l)m/j,„+i/j+ec,m-l [ «j+e,,m 
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Solving these equations, one obtains the frequencies u and the wave functions (wj.n, "^j.^) 
of excitations, which allow us to discriminate stability of a steady state 36|. When 
all the excitations satisfy the condition A/'w > 0, the steady state is stable, where 
A/" = Xli Xlnd'^^j:"!^ ~ I'^j.raP) is ^^^ normalization constant. If there exist excitations sat- 
isfying JVu < 0, the state is energetically unstable. This instability is called the Landau 
instability. The emergence of excitations with complex frequencies, i.e. |lm[a;]| > 0, sig- 
nals the DI which means exponential growth of the fluctuation in time. In a system of 
atomic gases at sufficiently low temperatures, while the LI cannot destabilize the system 



because of the 
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ack of energy-dissipation processes, the DI can drastically break down the 
38J . Hence, it is important to calculate the critical momentum for the DI, 



even though it is always larger than the Landau critical momentum. 



D. Local momentum 

We characterize the transport in trapped systems by the local momentum rather than 
the center-of-mass (COM) momentum, and compare it with the critical momenta in homo- 
geneous systems. The local momentum in the direction a is expressed as 



p? = sin-i I ^ X^^ I ' (9) 



I^n'in^,, 



where n? = |$jP is the condensate density at the site j and L" = ($*$j+eQ — '^j'^l+ea)/^ 
is the current which flows from the site j to the site j + e^. We will show in Sec. IIVI that 
significant damping of dipole oscillations occurs when the maximum local momentum in a 
trapped system exceeds the critical momentum for the DI in a homogeneous system. 

III. CRITICAL MOMENTA IN HOMOGENEOUS OPTICAL LATTICES 

In this section, we numerically solve Eqs. ([7]) and (|H]) to obtain the excitation spectra for 
the homogeneous Bose-Hubbard systems at i^ = 1 in ID, 2D, and 3D, where v is the filling 
factor. Using the discriminant described in the previous section, we determine the critical 
momenta for the LI and the DI. Notice that the critical momenta for the DI have been 
obtained in Refs. 16|, ll7| by solving the equation of motion (|3]) in real-time in a situation 
where the on-site interaction U or the flow momentum is slowly increased. A clear advantage 
of our method over the previous work is that the critical momenta for the LI is also available. 
We also emphasize that one can calculate the critical momenta for the DI more accurately. 
In the limit of t/ — )■ 0, for instance, while the critical momentum obtained in Refs. Il6 



is a little larger than the prediction by the GP mean-fleld theory, i.e. p = 7r/2 [39|, our 
method gives a precise agreement with the GP prediction as we will see below. 

To calculate the excitation spectra for moving Bose gases, we introduce the phase twist 
term in Eq. ([T]) as 

«j«j+e. -> ajaj+e.e^P'^'^ (10) 

which represents the situation where the optical lattice is moving at a constant momentum 
—p. Forcing the coefficients /j,„ to be real numbers, we carry out imaginary-time evolution 
of Eq. (j3]) to obtain the steady state with zero momentum, which is equivalent to a state 



with momentum p in a static optical lattice. Substituting the coefficients /j,„ of the steady 
state into Eqs. ([7j) and ([8]) and solving them, we obtain the excitation spectra of the current- 
carrying state. 
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ka/ 7T 

(d) 

FIG. 1: Excitation spectra fvjj{k) in ID homogeneous optical lattices for U/J = 3, ;/ = 1, and 
different values of p. (a) : p = 0, (b) : p = pi,i ~ 0.85, where the LI sets in, (c) and (d) : 
p = 0.94 >pDii. 

In Fig. [H we show the excitation spectra of Bose gases in ID homogeneous optical lattices 
for f/ = 3 and different values of the flow momentum p. We plot the ffist branch of the 
excitations that corresponds to the well-known Bogoliubov spectrum, because only this 
branch is relevant to the LI and the DI in most of the parameter regions. When p = 0, 
this branch has a phonon dispersion at |A;| ^ 1 (Fig. [U^a)). When p increases, the slope of 
the phonon dispersion for /c < decreases, and it reaches zero at p = pu — 0.85 as seen in 
Fig. [U^b). This signals the onset of the LI. With p increases further and exceeds a certain 
threshold, the excitations at \k\ <^ 1 have a finite imaginary part, which signals the DI. We 
call this threshold poii- 

In Fig. O we show the critical momenta for the LI (dashed red lines) and the DI caused 
by excitations with \k\ <^ 1 (solid black lines) in ID, 2D, and 3D. In any dimensions. 
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FIG. 2: (Color online) Critical momenta versus U/J in homogeneous optical lattices for v = 1. (a): 
ID, (b): 2D, and (c) 3D. The dashed red, solid black, and dotted blue lines represent the critical 
momenta for the LI (pLi)) the DI caused by excitations with long wavelengths (j'dii); and the DI 
caused by excitations with short wavelengths (pdk)- 



Pli — ^ and pdii ~^ 7r/2 when t/ — ?> 0. When U increases, pdii monotonically decreases 
and Pli approaches poii- In the strongly interacting region {U > f/c/2), Pli takes almost 
the same value as poii, and both momenta reaches zero at [/ = Uc-, where f/c denotes the 
SF-MI transition point. The behaviors of pdii are consistent with the previous work of 



Refs. 
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When the interaction is sufficiently weak in ID, a DI caused by excitations with short 
wavelengths precedes the DI caused by excitations with \k\ <^ 1. The critical momentum for 
this DI pdi2 is plotted by the dotted blue line in Fig. [2](a). To understand this instability, we 
show in Fig. [3] the excitation spectra including the Bogoliubov mode, the amplitude mode, 
and the anti-modes of the two modes 40|] for U = 0.1 and different values of p. When p 
increases, the amplitude mode declines and exhibits an avoided crossing with the Bogoliubov 
mode as seen in Figs. [3]^b) and (c). With increasing p further, the amplitude (Bogoliubov) 
mode is coupled with the anti-Bogoliubov (anti-amplitude) mode and these coupled modes 
cause the DI 36|. We suggest that this DI leads to the formation of density waves with 
ordering vectors k ~ iO.Svr and iO.Svr. While this type of DI has been previously found in 
the presence of off-site interactions 28|, l4l|-|43l|. this is the first example of such a DI in the 



standard Bose-Hubbard model only with the on-site interaction. 



-05 0.5 

ka/ 7T 

(a) 



S 

3 

£, -1 
<u 
a: -2 

-3 



-05 05 

ka/ 7T 

(b) 



0.05 
— , 0.04 
5 0.03 
J 0.02 
£ 0.01 




-0.5 0.5 

ka/ 7T 

(c) 



-0.5 0.5 

ka/ 7T 

(d) 



-05 05 

ka/ 7T 

(e) 



FIG. 3: Excitation spectra fujj{k) in ID homogeneous optical lattices for the weak interaction 
U/J = 0.1, V = 1, and different values of p. (a): p = 0, (b): p = 0.5, where the avoided crossing 
of the Bogoliubov and amplitude modes starts to occur, (c): p = 1.2, (d) and (e): p = 1.3, where 
excitations with short wavelengths causes the DI. The solid black, red dashed, and blue dotted 
lines represent the normal modes, the anti-modes, and the modes whose imaginary part is finite. 

IV. DIPOLE OSCILLATIONS IN THE STANDARD BOSE-HUBBARD SYSTEM 

In this section, we consider systems of Bose gases in ID and 2D optical lattices combined 
with a parabolic trap. We fix the total number of particles to be A^ = 45 in ID and A^ = 3000 
in 2D such that the Mott plateau at unit filling forms in the central region of the trap when 
U exceeds Uc- We first calculate the ground state via imaginary-time evolution of Eq. (J3]). 
We next calculate real-time dynamics subjected to a sudden displacement of the trap center 
that induces a dipole oscillation. Notice that while similar dynamics have been studied in 
Refs. 



35 



44l | using the same Gutzwiller approximation, the previous studies did not address 



the relation between the dipole oscillations and the critical momenta, which is our main 
interest here. 
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FIG. 4: Time evolution of the COM Xcom{t) in ID optical lattices combined with the paraboric 
trap for U/J = 6. (a): stable dipole oscillation at D = 4.3, (b): damped dipole oscillation at 
D = 4.A. 

A. Stability and damping of dipole oscillations 

To illustrate basic properties of the dynamics of dipole oscillations, we start our analyses 
with the ID case. In Fig. HI we show typical examples of the COM motion for stable 
and damped dipole oscillations, where the COM is given by Xcom = J2jJ^j/ J2j^j- The 
amplitude of the dipole oscillation hardly changes as long as the displacement D is smaller 
than a certain critical value Dc as seen in Fig. ID^a). This dissipationless motion is a clear 
characteristic of superflow. In contrast, the oscillation is remarkably damped when D^ is 
exceeded (see Fig. |4](b)). It is worth noting that the dynamics for these two cases can be 
clearly distinguished because the transition to the damped oscillation occurs very sharply. 
Thanks to this property, one can determine the critical displacement accurately. 

B. Local momentum 

For comparison with the critical momenta in the homogeneous system, we calculate 
two types of momentum in the trapped system, namely the COM momentum, Pcom = 
sin~^(t;com/2), and the local momentum defined in Eq. (jH]), where t^com is the COM velocity. 
In Fig. [5], we show the snap shots of the local momentum pj, the condensate density n^ = 
|$jp, and the density rij at the critical displacement D = D^ for a weak interaction U = 2 
(a) and a strong interaction U = 10 (b). The time is set such that the momentum of COM 
takes the maximum value during the time evolution. We define this time as tmax- For the 
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FIG. 5: (Color online) Local momentum pj (solid black line), local condensate density n^ = |$jp 
(dashed red line), local density rij (dotted black line), local current /_,• (dot-dashed green line) in 
ID at D = Dc and t = i^ax, where U/J = 2 (a) and U/J = 10 (b). 



weak interaction, the local momentum is almost constant. On the other hand, for the strong 
interaction in Fig. Mjo), the local momentum is significantly dependent on the position so 
that it peaks at the unit-filling points. This happens because the condensate density in 
the regions close to unit filling is strongly suppressed as a precursor of the formation of the 
Mott plateaus. Due to the strong spatial dependence of the local momentum, the COM 
momentum is noticeably smaller than the maximum local momentum, especially near the 
Mott transition. 

C. Critical momentum 

In Fig. [6l we show the maximum local momentum Pmax and the COM momentum Pcom at 
D = Dc and t = tmax as functions of U. For comparison, we also plot the critical momenta 
Pdii obtained in the homogeneous systems at unit filling. We find that Pmax(^ = ^max) in 
the trapped system quantitatively agrees with poii- In contrast, Pcom(t = tmax) drastically 
deviates from poii- These results lead to the important conclusion; for detecting the critical 
momenta for the homogeneous lattice systems from the dipole oscillations, one needs to 
measure the local momentum rather than the COM momentum that can be measured rela- 
tively easily from the time-of-fiight images 
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13| . While the local momentum has never 



been experimentally measured so far, it should be available in future experiments because 
recent experiments have developed techniques to address local observables at the scale of a 
single lattice site [TI. l45l| . 
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FIG. 6: (Color online) Maximum local momentum Pmax (solid black lines with symbols) and the 
COM momentum Pcom (dashed red lines) at D = Dc and t = tmax as functions of U/J. (a) : ID 
and (b) : 2D. The dotted blue lines represent pBii in the homogeneous lattice systems, which are 
plotted in Figs. [21^a) and (b) with the solid black lines. 

D. Mode coupling 
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FIG. 7: The motion of the COM a;com(0 with damping and revival in ID optical lattices combined 
with a parabolic trap for U/J = 7.78. 

We now focus on long-time dynamics of the dipole oscillation. In this case, we find that 
the dipole oscillation with the damping and the revival occurs below the critical momentum 
in a particular range at the interaction U = 7.5 — 7.8. Similar behaviors have been found in 
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FIG. 8: (Color online) (a) Frequencies of the first mode oji (solid black line) and the second mode 
^2 (dashed red line) in ID lattices combined with a parabolic trap, (b) Ratio 002/101. 



Refs. (35|, |4J|. In Fig. [TJ we show this unusual oscillation for U = 7.78. For understanding 
this phenomenon, we calculate the excitation spectrum by diagonalizing Eqs. ([7]) and ([8]). 
In Fig. [8]^a), we plot the frequencies of the first and second modes, which correspond to the 
dipole mode and the breathing mode, respectively. We also show the ratio of the breathing- 
mode frequency U2 to the dipole-mode frequency Ui in Fig.[H](b). We find that the ratio uj^I'^x 
is nearly equal to two in the regime where the damping and the revival occur. Therefore, this 
phenomenon can be interpreted as indicating that dipole and breathing modes are coupled 
by the nonlinear effect in Eq. ([3]). In this regime, although the amplitude of the dipole 
oscillation decreases due to this mode coupling even below the critical momentum, we can 
determine the critical momentum because the damping caused by mode coupling is much 
smaller than that caused by the DI. 



V. DIPOLAR HARDCORE BOSONS 



In this section, we consider dipolar bosons confined in 2D optical lattices combined with 
a parabolic trapping potential. For simplicity, we consider the hardcore limit (f/ — )■ 00), 
in which the local Hilbert space is spanned only with |0) and |1), and assume that the 
dipoles are polarized to the direction perpendicular to the lattice plane. This system is well 
described by the following Bose-Hubbard model. 
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H = -jJ2J2(^lai+-^ + h.c.) + -Y, ^3 
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where V is the strength of dipole-dipole interaction, rj i = a|j — 1| is the distance between 



the sites j and 1. Since the dipolar interaction potential decays as 



r-~j 7" 



j,l ' 



we only include 



the dipolar interactions in the range rj i < 7a. Introducing this cutoff merely changes some 
quantitative natures of the system very slightly, but does not affect any qualitative natures, 
such as the quantum phase diagram and the critical momentum, when V^ < 10. We analyze 
this system by using the Gutzwiller approximation as in the previous section. We again 
adopt the units oih=J = a = l. 

A. Ground state 
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FIG. 9: (Color online) The average density rictr (dotted blue line), The average condensate density 



n" 



(solid black line), and the static structure factor S'ctr('?r, vr) devided by 5 (dashed red line) 



around the trap center in the ground state of Eq. (fTTj) are plotted as functions of V/J, where 
n = O.OIJ and N = 900. 

Let us calculate the ground state of Eq. flTTl) . We fix the curvature of the parabolic 
potential and the total number of particles to be 1) = 0.01 and A^ = 900. In this case, the 
ground state can be the following two types of state as long as V^ < 4.5. For small V, the 
dipolar gas in the ground state consists only of a SF region. When V is sufficiently large, 
the ground state has a region of the SS with the checkerboard density wave order around 
the trap center, and this SS region is surrounded by a SF region. As a metaphor, one may 
imagine a sunny-side up whose yolk and white correspond to the SS and SF regions. In order 
to identify the transition between the two states, we show in Fig. [9] the average condensate 
density and the static structure factor around the center of the trap, which are defined as 
Eiectr l-^jlVSS and S,U^,7t) = Enectr ^^^"'"^■^J-'^K^;)/25. Here ^^^^^ denotes the 



n 



ctr 



summation of the 5x5 sites around the center of the trap. For comparison with the results 
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in Ref. [28|, we also plot the average density rictv = (^o + neJ/2. By using these quantities, 



we identify the state with ra^^-j, ^ and Sdri'^,'^) = as the SF state and the state with 
n^tj. 7^ and S'ctr(7r,7r) 7^ as the SS state. We find that the boundary of the SF and the 
SS is 1/ = 3.26. 

B. Critical momentum 

Having obtained the ground state of Eq. (11 II) . we next calculate the dipole-oscillation 
dynamics following the way illustrated in the previous section. We find that the critical 
displacement Dc can be accurately determined also in the case of dipolar hardcore bosons. 
By the solid black line with symbols in Fig. [TOT a). we show the maximum local momentum 
Pmax ai D = Dc and t = tmax as a function of V. We clearly see that pmax agrees well with 
the critical momentum for the DI in the homogeneous lattice system at u = OA obtained in 



Ref. [28|]. Pmax takes a distinct minimum at a point close to the transition to the SS state 
and is significantly smaller in the SS state than in the SF state. We emphasize that the 
agreement with the homogeneous case is found in the SS state even though there also exists 
the SF region with a different critical momentum. This can be attributed to the fact that 
the critical momentum in the SS phase is smaller than that in the SF phase and that the 
maximum local momentum is always taken at the center of the SS region where the local 
condensate density is smallest. 

Let us discuss the feasibility of identifying the SS state from the critical momentum. 
If one can measure the local momentum, the above-mentioned properties of the critical 
momentum in the SS state allow for distinguishing the SS state from other possible states, 
such as the SF, MI, and density-wave insulating states. However, since it should be easier 
to measure the critical displacement rather than the critical local momentum, we show the 
critical displacement as a function of V in Fig. [TUT b). where the significant reduction of D^ 
in the SS state is clearly seen. Thus, we suggest that it is possible to identify the SS state 
by measuring the critical displacement Dc. 
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FIG. 10: (a) The solid line with symbols represents the maximum local momentum Pmax at D = Dc 
and t = imax hi the system of Eq. ([TT]) as a function of V/J. The dashed red line represents the 



critical momentum for the DI in the homogeneous case at z/ = 0.4 
Dc as a function of V/J. 



28[. (b) Critical displacement 
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FIG. 11: (Color online) Snap shots of the density distribution rij during the dipole-oscillation 
dynamics of Eq. ^ with V/J = 3.2, where Jt/h = (a) and Jt/h = 310 (b). 



C. Dynamical transition 
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that when the superflow 



It has been predicted for the homogeneous case in Ref. 
momentum increases in the SF state near the phase boundary to the SS state, the transition 
to the SS phase can occur due to the DI caused by excitations with k = (7r,7r). Since 
typical experimental setup includes a parabolic trapping potential, it is important to confirm 
whether this dynamical transition induced by superflow occurs or not in trapped systems. 
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FIG. 12: Time evolution of the structure factor ^'(vr, vr) during the dipole oscillation, where V j J = 
3.2 and D = 2. 

For this purpose, we choose V = 3.2 such that the ground state is in the SF state close to 
the transition point. In Fig. [TTT a). we show the density profile of the ground state. Taking 
this state as the initial state, we compute the dipole-oscillation dynamics of Eq. (ITT]) setting 
D = 2. To quantify the formation of the density wave order associated with the dynamical 
transition to the SS, we compute the time evolution of the static structure factor 5'ctr(7r, vr), 
which is shown in Fig. [121 Sctri'^i'^) exhibits repetitive growth and collapse after the time 
t ^ 150. In Fig. [TlTb). we show the density distribution at t = 310, where the growth of the 
S'ctr(7r, vr) is the most prominent within the time scale plotted in Fig. [121 There we clearly 
see that the checkerboard density wave forms around the trap center. Thus, the dynamical 
transition to the SS can occur in the dipole-oscillation dynamics in the trapped system. 

VI. SUMMARY 

We have studied the superfluid transport of Bose gases in optical lattices by using the 
Gutzwiller approximation. In one-dimensional (ID), 2D, and 3D homogeneous systems, we 
determined the critical momenta for the Landau instability (LI) and the dynamical instabil- 
ity (DI) from the excitation spectra. Especially, we have found a DI caused by excitations 
with short wavelengths in ID when the on-site interaction is very small. In trapped system, 
we have analyzed the dynamics of dipole oscillations induced by suddenly displacing the trap 
center. We have found that the critical momentum defined by the maximum local momen- 
tum in the trap system quantitatively agrees with that in the homogeneous system at unit 
filling. We have also shown that a resonance phenomenon is caused by the coupling between 
the dipole-oscillation mode and the breathing mode. Moreover, we have investigated the 
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critical momentum of dipolar hardcore bosons and shown that the critical momentum of the 
SS state is smaller than that of the SF state as in the case of the homogeneous system. This 
result allow us to suggest that the SS state should be identified by measuring the critical 
displacement. We finally showed that the dynamical transition from SF state to SS state 
can be induced by dipole oscillations. 
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